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Q . Deep long-slit spectra of the diffuse ionized gas halos of the edge-on spiral 

I . galaxies NGC 4302 and UGC 10288 are presented. Emission lines are detected 

up to about z = 2 kpc in NGC 4302, and to nearly z = 3 kpc on the north side of 
^ . UGC 10288. For both galaxies, the line ratios [N Ilj/Ha and [S Il]/Ha increase 

^ ! with z in accordance with dilute photo- ionization models. Runs of [S 11] /[N II], 

if^ . and for UGC 10288, the run of [O III] /Ha, however, are not explained by the 

^ I . models. Scale height determinations of their DIG halos are generally lower than 

O ■ those of galaxies with more prominent extraplanar DIG features. 

^ , These data, along with previously presented data for NGC 5775 and NGC 

P_|. 891, are used to address the issue of how DIG halos are energized. Composite 

Q . photo-ionization/shock models are generally better at explaining runs of line 

^ , ratios in these galaxies than photoionization models alone. Models of line ratios 

d '. in NGC 5775 require a greater contribution from shocks for filamentary regions 

J> ! than for non-filamentary regions to explain the run of [O III] /Ha. In either case, 

^ ! the [S II] /[N II] ratio is not well fit by the models. Composite models for UGC 

10288 are successful at reproducing the run of [S II] /[N II] for all but the the 
highest values of [N II]/Ha; however, the run of [O IIl]/Ha vs. [N II]/Ha does 
not show any discernible trend, making it difficult to determine whether or not 
shocks contribute to the layer's maintenance. 

We also examine whether the data can be explained simply by an increase in 
temperature with 2; in a pure photo-ionization model without a secondary source 
of ionization. Runs of [S Il]/Ha, [N II]/Ha, and [S II]/[N II] in each of the four 
galaxies are consistent with such an increase. However, the runs of [O III] /Ha 
vs. z in NGC 5775 and UGC 10288 require unusually high ionization fractions 
of O^^ that can not be explained without invoking a secondary ionization 



^Visiting Astronomer, Kitt Peak National Observatory, National Optical Astronomy Observatories, 
Tucson, AZ 



- 2 - 



source or at the very least a much higher temperature for the [O III] -emitting 
component than for the [S II]- and [N II] -emitting component. An increase in 
temperature with z is generally more successful at explaining the [O III] /Ha run 
in NGC 891, with the ionization fraction of O"'""'" remaining relatively low and 
nearly constant with z. 

Subject headings: diffuse radiation — galaxies: ISM — galaxies: spiral — stars: 
formation 



1. Introduction 

A vertically extended layer known as the Reynolds Layer or the Warm Ionized Medium 
(WIM) of the Milky Way contains the majority of the ionized gas in the interstellar medium. 
The layer fills about 20% of the ISM volume with a local midplane density of ~0.1 cm~^ 
and a scale height of ~1 kpc (Reynolds 1993). Energetic arguments favor photo-ionization 

by massive stars in the disk for the layer's maintenance, as roughly 15% their ionizing 
output is sufficient to maintain the layer in the solar neighborhood (Reynolds 1993). The 
Wisconsin Ha Mapper (WHAM; Reynolds et. al. 1998b) has been instrumental in recent 
years in characterizing the WIM. In external galaxies, where the layer is more commonly 
referred to as Diffuse Ionized Gas (DIG), analogous information has been obtained through 
the use of narrow-band emission-line imaging, long-slit spectroscopy, and Fabry-Perot 
observations [e.g. Rand, Kulkarni, & Hester 1990 (hereafter RKH); Rand 1996, 1997, 1998; 
GoUa, Dettmar, & Domgorgen 1996; Ferguson, Wyse, & Gallagher 1996; Wang, Heckman, 
& Lehnert 1997; Hoopcs, Walterbos, & Rand 1999; CoUins et. al. 2000 (hereafter CRDW)]. 
Estimates of the diffuse fraction in DIG layers put the ionization requirement closer to 40% 
of the massive stars' ionizing output (e.g. Hoopes et. al. 1996; Ferguson et. al. 1996). The 
advantage in observing external edge-on galaxies is the ability to obtain an overview of the 
entire halo, so that brightness, spatial distribution, and emission line ratios can be better 
characterized as they vary with z. 

One of the most important issues regarding these layers is their ionization. In the 
WIM, as well as in external DIG layers, [S II] A6716/Ha and [N II] A6583/Ha are enhanced 
relative to [H II] regions [e.g. Haffner, Reynolds, & Tufte 1999 (hereafter HRT); Ferguson, 
Wyse, & Gallagher 1996; Rand 1997]. Such a scenario is predicted by photo-ionization 
models in which the dilution of the radiation field, measured by the ionization parameter, 
[/, increases as photons propagate away from H II regions [Domgorgen & Mathis 1994; 
Sokolowski 1994 (hereafter S94); Bland-Hawthorn et. al. 1997]. Species such as S which are 
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predominantly doubly ionized in H II regions are able to recombine to a singly ionized state 
as the field becomes increasingly dilute. Such an effect shoTild not be as significant for N or 
O, which have considerably higher second ionization potentials and thus may not be fully 
doubly ionized in H II regions. 

The photo-ionization scenario is not without problems, however. Firstly, in an 
increasingly dilute field, photoionization models predict that the ratio [S II] /[N II] should 
increase with z due to differing second ionization potentials of S and N. However, [S II] /[N II] 
vs. z has been found to remain nearly constant in the WIM (HRT) and NGC 891 (Rand 
1998), and in the case of NGC 4631, the variations with z are not as pronounced as those 
predicted by the models (GoUa et al. 1996). In addition, the ratio [O III] A5007/Hq! 
should decrease with z in a photo-ionization scenario as the field becomes increasingly 
dilute. The opposite trend, however, is seen in NGC 891 (Rand 1998), where [O III] /Ha 
increases significantly from the midplane to 2; = 2 kpc, as well as in NGC 5775 (Rand 2000; 
Tiillmann et. al. 2000). [O III] /Ha enhancements in DIG relative to H II regions are also 
seen in a number of more face-on galaxies, including M51 (Wang et. al. 1997). Finally, 
photoionization models have trouble explaining values of [N II] /Ha and [S II] /Ha that are 
often 1 (Wang et. al. 1997; Rand 1998; HRT). Models of S94, who took observed line 
ratios in NGC 891 into consideration, can only explain such ratios by assuming an upper 
mass cutoff for the stellar IMF of 120Mq, hardening of the radiation field as it propagates, 
and a depletion of important gas-phase coolants onto dust grains. Observations of the weak 
He I A5876 hue, however, in the Reynolds Layer (Reynolds & Tufte 1995) and in NGC 891 
(Rand 1997) each imply a significantly softer spectrum than indicated by the forbidden 
hues. 

One possibility that has been forwarded to explain some of this line ratio behavior is 
to assume a secondary ionizing source that contributes a small but increasing fraction of 
the Ha emission with z. Such secondary sources include shocks (ShuU & McKee 1979) 
or turbulent mixing layers (TMLs; Slavin, ShuU, & Begelman 1993), each of which may 
produce large amounts of [O III] emission. Other, less well-studied ionization sources may 
be relevant in halos, such as ionization from soft X-rays or cosmic rays. 

Another approach that has been taken is to assume, instead of the ionization 
parameter or state changing with ^, that the line ratio behavior can be explained by the 
gas temperature increasing with z. This requires a source of additional heating, such as 
photo-electric heating from dust grains (Reynolds & Cox 1992) or dissipation of interstellar 
turbulence (Minter & Balser 1997), that becomes more significant at lower gas densities, 
but may not require any additional source of ionization. A temperature increase with z 
is successful in explaining enhancements in [N II] /Ha and [S Il]/Ha with z, as well the 
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constancy in [S n]/[N II] in the Reynolds Layer up to z = 1.75 kpc (HRT; Reynolds, 
Haffner, & Tufte 1999). Such a rise in temperature could also explain an increase in 
[O III] /Ho;, though data on the run of [O III] /Ha vs. z in the WIM are currently lacking. 

The goal of this research is to attempt to better understand how DIG halos are 
energized. We attempt to determine whether measured runs of line ratios versus z can be 
explained by: (a) photo- ionization with a contribution from shock ionization that increases 
with height above the midplane, and/or (&) only photo-ionization, but with an increase in 
gas temperature with z. We examine our data in the context of each of these scenarios, 
realizing that elements of both could be important for understanding the ionization of DIG 
halos in general. To this end, we have obtained data for two galaxies: NGC 4302 and 
UGC 10288. We also employ data for NGC 5775, already presented by Rand (2000). For 
the purpose of assessing the rising temperature scenario, we use data for NGC 891, which 
has been compared with composite shock/photo-ionization models previously (Rand 1998). 
The observations of these galaxies are discussed in § 2. Line ratio data are presented in 
§ 3 along with comparisons to composite shock/photo-ionization and rising temperature 
models. We briefly summarize our results in § 4. 

2. Observations 

Observational details for the NGC 891 (Rand 1998) and NGC 5775 (Rand 2000) data 
have been discussed in previous papers, though some modifications have been made which 
will be discussed in this section. Spectra for UGC 10288 and NGC 4302 were obtained 
at the KPNO 4-m telescope on 1999 June 10-13 and 1998 February 23-26, respectively. 
Slits for all four galaxies run perpendicular to the plane. The slit position for NGC 4302, 
covering the most prominent region of extraplanar DIG, is shown in Figure 1, overlaid on 
the Ha image. The Ha image of UGC 10288, showing httle extraplanar Ha emission, is 
shown in Figure 2 with slit position overlays. The slit positions for UGC 10288 were chosen 
to be roughly 30" east and west of the galaxy's center. Due to the presence of a foreground 
star on the axis perpendicular the galactic disk 30" west of the radio continuum center, slit 
2 was placed at a position 35" west of center. A summary of the observations is given in 
Table 1. 

The observations of NGC 4302 utilized the KPC-24 grating in conjunction with the 
T2KB 2048x2048 CCD, providing a dispersion of 0.5lA per pixel, a resolution of 1.3A and 
a useful coverage of about 800A. The slit length is 5', the slit width was set at 3", and the 
spatial scale is 0.69" per pixel. For red line observations, the grating was tilted to give an 
approximate central wavelength of 66OOA, allowing Ha, [N II] AA6548, 6583, [S II] AA6716, 
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Fig. 1. — The slit position for the observations of NGC 4302 overlayed on the Ha image of 
Rand (1996). The image does not show the entire 5' sht length. 
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6731, and [O I] A6300 to be observed. Due to poor weather during the 1998 observing run, 
spectra containing [O III] A5007 could not be obtained. 

For observations of UGC 10288, the KPC-007 grating was used with the T2KB 
2048x2048 CCD, providing a dispersion of 1.42A per pixel, a resolution of S.sA and a useful 
coverage of about 2100A encompassing the spectral lines Ha, [N II] AA6548, 6583, [N II] 
A5755, [S II] AA6716, 6731, He I A5876, [O I] A6300, [O III] A5007, and H/3. The slit length 
is 5', the slit width was set at 2", and the spatial scale is 0.69" per pixel. 

Data reductions were performed with the IRAFQ package. Projector flats were used to 
remove small-scale variations in response, while sky fiats were used to determine the slit 
illumination correction. Arc lamp exposures were used to calibrate the wavelength scale 
as a function of location along the slit, and the spectral response function was determined 
with standard stars. Regions of the spectra free of emission from the galaxies were used 
to form sky-subtracted spectra. In order to increase sensitivity for the determination of 
line parameters, blocks of 5 pixels were averaged in the spatial direction. Since Gaussians 
appear to accurately represent the line profiles, Gaussian fits were used to determine line 
properties. Uncertainties in line parameters for NGC 4302 were determined with the IRAF 
task splot. The use of this method for NGC 5775 and UGC 10288 yielded error bars that 
were far greater than uncertainties suggested by the data. In this case error bars were 
estimated by determining fiux variations along the slit for sky lines of a range of intensities, 
from which a relation between these variances and sky line intensity could be established. 
This relation was then used to determine fiux uncertainties for the various emission lines. 
Uncertainties for Ha line strengths in NGC 5775 also take into account the coincidence 
of the Ha line with a blended pair of sky lines. This factor was not considered in the 
presentation of NGC 5775 line data in Rand (2000). 

In order to compare [O III] A5007 emission to that of Ha it was necessary to make 
a reddening correction to the [O III] line in NGC 5775, UGC 10288, and NGC 891. This 
correction was not performed for the presentation of the [O III] /Ha ratio in Rand (2000) 
for NGC 5775. In all three galaxies, H/5 could be measured for all points where [OIII]A5007 
could be measured, thus allowing extinction corrections to be made for the full range in 
z. The correction was determined by assuming a uniform foreground dust screen and 
comparing measured values of the H/5/Ha ratio to the value calculated by Osterbrock 
(1989) for a 10^ K gas. The assumption of a 10^ K gas is not contradictory to a scenario 
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including a rising temperature with z in that the H/3/Hq; ratio is not nearly as temperature 
sensitive as the ratios of forbidden hncs to Balmer hnes: from 5,000 to 10,000 K, H/3/Hq; 
increases by 5%, whereas the ratios of [N II], [S II], and [O III] to Ha increase by an order 
of magnitude. In NGC 5775, the correction to the [OIII] flux decreases from a factor of 2.9 
in the midplane to no measurable extinction by z = 3 kpc. The correction is more variable 
for UGC 10288, with a maximum correction of 3.6 at ^ = 500 pc on the south side. In 
NGC 891, the correction decreases from 2.2 in the midplane to 1.5 at z = 1.6 kpc. In both 
NGC 5775 and NGC 891, the H/3 line is superposed on a weak stellar absorption line close 
to the midplane. However, the absorption line is much broader than the emission line and 
thus can be accounted for by the background fitting process in the line fitting routine. 

3. Results and Discussion 

3.1. Line Data and Line Ratios 

The Ha, [N II] A6583, and [S II] A6716 lines are detected up to nearly z — 2 kpc on 
either side of the plane of NGC 4302. The [S II] A6716 line was not detected &t z — 1.6 
kpc on the west side, though the fainter line at 6731 A was detected, due to confusion 
with noise in that part of the spectra. The flux of the 6716A line at that part of the 
spectra was inferred from that of the 6731 A line assuming the low-density limit of [S II] 
A6716/A6731 = 1.4. The vertical runs of [N II] A6583/Ha, [S II] A6716/Ha, and [S II] 
A6716/[N II] A6583 are shown in Figure 3. The data are averaged over 5 spatial pixels. 
The [N Il]/Ha and [S ll]/Ha each show a similar rise with with the [N ll]/Ha rising 
from about 0.4 in the midplane to 1.0 at 2; = 1.2 kpc on the east side, and to nearly 1.4 at 
z = 1.5 kpc on the west side. The [S n]/[N II] remains fairly constant at about 0.6 for the 
full range of z, with slight increases in the midplane and beyond z — 1 kpc. 

In UGC 10288, the Ha, [N II] A6583, and [S II] A6716 lines are all detected well above 
the midplane, with the [N II] A6583 and Ha lines being detected up to nearly z = 3 kpc on 
the north side for both slit positions. Detections of the [O III] A5007 and [O I] A6300 lines 
are limited to z <^ 1 kpc, except for the south side of slit 1 where both lines are detected up 
to z = 1.5 kpc. The vertical runs of [N II] A6583/Ha, [S II] A6716/Ha, [S II] A6716/[N II] 
A6583, [O III] A5007/Ha, and [O I] A6300/Ha for both slits are shown in Figure 4. The 
data are averaged over 5 spatial pixels. For sht 1, [N II]/Ha rises from 0.3 at the midplane 
to 1.6 at 2; = 2 kpc on the north side and 1.3 at 2; = 2 kpc on the south, while for slit 2 the 
ratio remains somewhat lower, rising to 1.2 for z — 1.5 kpc on either side of the disk. The 
[S n]/[N II] ratio shows a similar range of values for each slit, varying from 0.8 to 1.4 with 
the lower values generally occurring in the midplane. These values are significantly higher 



-8- 



than those of the other three galaxies. The [O III]/Hq; ratio shows some evidence for an 
increase with z on the south side of either sht, reaching values as high as 0.75 in Slit 1 and 
0.4 in Slit 2. On the north side of each slit, values reach no higher than 0.5 in Slit 1 and 
0.4 in Sht 2. Such values of [O III]/Hq; are consistent with those observed at comparable 
^-heights in NGC 5775. 

Line ratio data for NGC 5775, except for the [O I] A6300/Hq; and extinction corrected 
[O III] A5007/Ha ratios, have been presented previously (Rand 2000). These, along with 
[N II] A6583/Ha, [S II] A6716/Ha, and [S II] A6716/[N II] A6583 ratios are shown in Figure 
12. The data are averaged over 10 spatial pixels. An increase with z in the extinction 
corrected [O III] /Ha for the NE side of Slit 1 and the SW side of Slit 2 is still indicated 
by the data, where values range from 0.2-0.4 in the disk to 1.1 in the halo. The [O I/Ha 
ratio ranges from about 0.02 in the midplane to 0.06 z — 2.5 kpc in Slit 2, while the 
[O I] A6300 line was not detected for Slit 1. Because of the 85° (Irwin 1994) inclination of 
the galaxy, points within 10" of the midplane {z ^ 1200 pc in the figures) on the spatial 
axis reflect in-plane, highly inclined disk structure rather than true vertical structure. In 
these lines of sight, H II regions and areas between them have a greater effect on line ratio 
variations than extraplanar diffuse gas. This fact explains why some line ratio minima are 
not at z = kpc. 

Line ratio data for NGC 891, previously presented by Rand (1998), can be seen in 
Figure 13. 

3.2. DIG scale heights 

In this section we attempt to determine the Ha emission scale height. Hem, of the 
galaxies' DIG halos at the various slit locations. We attempt to model only an exponential 
halo component, avoiding any disk component where contamination by bright H II regions 
could complicate the analysis. We plot the logarithm of Ha intensity versus z for each 
slit, then fit an exponential component with a least-squares fit except in cases where an 
exponential is obviously a bad description. Plots of logarithmic Ha intensity vs. z, along 
with exponential fits, for NGC 5775, NGC 4302, and UGC 10288 are shown in Figures 5, 6, 
and 7, respectively. Results of the fitting procedure are summarized in Table 2. 

Position 1 of NGC 5775 covers the NE filament, DIG emission from which has been 
modeled previously (CRDW). The value Hem = 2300 pc is somewhat consistent with the 
value oi H — 3.0 kpc determined from image data. The image data however does include 
some [N II] contamination. The previously determined value assumes an [N Il]/Ha ratio 
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that does not vary with z. Since our data shows that this ratio in fact rises with z on the 
NE filament, from 0.5 at ^ = 2 kpc to 0.75 at z = 6 kpc, the value H — 3.0 kpc becomes an 
upper limit on the actual scale height. 

The DIG halos of NGC 4302 and UGC 10288 have been modeled previously as 
well (Rand 1996). The sht position for NGC 4302 coincides with one of the previously 

modeled regions, the scale height of which was determined to be = 525 pc (the image 
is not contaminated by [N II] emission). The values determined from spectroscopic data, 
Hera = 700 pc and Hem = 550 pc for the east and west sides, are close to those determined 
from image data. The scale height determined from the image, however, is an average 
over a 3 kpc radial extent of the disk and thus should not be in exact agreement with 
the spectroscopically determined value. Image data of UGC 10288 revealed very little 
extraplanar Ha emission, and thus profiles were not well modeled with exponentials. 
Increased sensitivity of the spectroscopic data have allowed the detection of Ha emission to 
greater z-heights. The resulting Ha profiles seem to be well modeled by exponentials with 
small scale heights on the south side, while more complicated structure is present on the 
north side. 

NGC 5775, having the most prominent extraplanar DIG emission, also has the DIG 
halo with the greatest scale height. In contrast, UGC 10288, showing very little extraplanar 
Ha emission, has a DIG halo of relatively small scale height. NGC 4302, having a number 
of extraplanar plumes of Ha emission, has a DIG scale height intermediate between that of 
NGC 5775 and UGC 10288. The halo of NGC 891 has been modeled previously (RKH), 
the scale height of which has been measured to be i?em = 2.4 kpc on the west side and 
Hem — 3.2 kpc on the east side from imaging data. These results suggest that not only do 
galaxies such as NGC 5775 and NGC 891 have brighter DIG halos, where sensitivity limits 
allow faint emission to be detected at much greater heights than for a galaxy with a lower 
Ha surface brightness, but the halos also have much larger scale heights. This adds further 
evidence to the notion of a dynamic halo where higher levels of star formation in the disk 
result in material driven further into the halo (e.g. Rand 1996; Hoopes et. al. 1999). 

3.3. Shocks as a Secondary Source of Ionization 

Some of the line ratio behavior observed in NGC 5775 and UGC 10288, namely the 
rise in [N II] /Ha and [S II] /Ha with z, is predicted by the photo-ionization models (e.g. 
S94; Bland-Hawthorn 1997). As previously stated, rises in these ratios are expected in such 
models as the radiation field becomes more dilute with z. If the field is in fact becoming 
more dilute (low values of the ionization parameter, U, which measures the "diluteness" of 
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the radiation field and is proportional to the ratio of ionizing photon number density to 
gas density) with then the models predict that the [O III] /Ha ratio should decrease 
away from the midplane. However, in NGC 5775, and to a lesser extent UGC 10288, we 
see a definite incease in [O III] /Ha with indicative of higher excitation conditions above 
the disk in such models (high-f/). Recent spectroscopy of NGC 5775 by Tiillmann et. 
al. (2000) also reveals a significant rise in [O III]/Ha with z, though these data are not 
extinction corrected. It is then possible that the [O III] emission originates from a different 
DIG component than that of the [N II] and [S II] . Consequently, a secondary mechanism, 
such as shock ionization, could account for the bright [O III] emission above the plane. 

In this section we explore the possibility that shock ionization is responsible for some 
of the DIG emission in these galaxies. One could also consider TMLs, since they also 
produce large amounts of [O III] emission. In fact in NGC 891, TMLs do just as well as 
shocks in fitting the observed data (Rand 1998). Our approach is that of Rand (1998), with 
some minor changes. We adopt the S94 matter-bounded photoionization model with the 
lowest terminal hydrogen column density, 1 x 10^^ cm~^, of individual gas clumps. The S94 
models feature a hard stellar spectrum (an upper mass cutoff of I2OM0), hardening of the 
radiation field as it propagates through the ISM, and depletion of gas-phase coolants onto 
dust grains. We do not consider the S94 radiation-bounded model as this model predicts 
extremely bright [O III] emission. Any composite shock model utilizing radiation-bounded 
models requires a very unusual secondary source of Ha emission. In the case of NGC 5775, 
such a scenario would require a secondary source that accounts for nearly half of the Ha 
emission at the midplane (assuming log U ~ —3.0 at 2; = kpc), and produces virtually 
no [O III] or [N II] emission. In addition, such a source must produce significant [O III] 
and [S II] emission at high- 2;, and still contribute at least 50% of the total Ha emission. 
Such emission characteristics could be explained with low-speed shocks near the midplane 
with a shock velocity increasing with z. However, the problem is further complicated by 
the fact that low-speed shocks produce copious amounts of [O I] emission relative to [N II]. 
In addition, basic energetic arguments render such a scenario, where the secondary source 
contributes half of the layer's ionization, highly unlikely as the ionizing flux from massive 
stars exceeds supernovae power output, presumably the initiator of these shocks, by a factor 
of 6 or 7 (Reynolds 1984, 1992). In any case, parameters describing these shocks require 
considerable tweaking to account for the observed run of line ratios. Matter-bounded 
models on the other hand, better match the data before composite modeling, and require a 
less complicated parameter space describing the secondary source. 

To obtain line ratios due to shock ionization, we use the models of Shull & McKee 
(1979; hereafter SM79). Line ratios in these models, especially [O IIl]/Ha, are highly 
dependent on shock velocity. Other variables affecting line ratios include the preshock gas 
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density and ionization state, abundances, and transverse magnetic field strength. These 
models assume the gas to be initially neutral at no = 10 cm~^, then penetrated by a 
precursor ionization front. More appropriate to the case of a partially photo-ionized DIG 
layer, would be a initial density of order hq ~ 0.1 cm^^, and a high initial ionization 
fraction. We mainly consider the models with shock velocities of Vs = 100 km s~^ and 

= 110 km in that these higher velocity models are best able to produce large amounts 
of [O III] emission, while keeping [O I] low relative to [N II]. SM79 also calculate a single 
model containing depleted abundances with Vg — 100 km s~^, which we also consider in the 
modeling to allow a comparison between shock and photo-ionization models with depleted 
abundances. In general, though, we do not consider this analysis as providing a definitive 
measure of the shock speed in these composite models given the simplicity of the model 
and the many unrealistic parameters which describe the shocks. 

For the composite modeling, we assume that some fraction of Ha emission arises from 
shock ionization, with that fraction possibly changing with z. We still assume a stellar 
radiation field distinguished by a decrease in ionization parameter, U, with z. To simplify 
matters we assume a single shock velocity for each composite model. We attempt to fit the 
composite model to the data by varying the percent contribution at a given value of U (and 
thus z). For NGC 5775 the modeling worked best by allowing shock contributions to begin 
at log U — —4.0, while models with shocks beginning at log U — —3.5 worked best for UGC 
10288. Such values are somewhat lower than would be expected for a \ow-z radiation field. 
However, kinematics of edge-on galaxies indicate that at low- 2;, the line emission we observe 
is mostly associated with the outer disk regions due to an absorbing dust layer. Since star 
formation is typically concentrated in the inner disk, the radiation field in the outer disk 
may be relatively dilute, in which case, the value of U at low- 2; may be quite low (Rand 
1998). 

We do not carry out a statistical test of the goodness of fit since our goal is to 
qualitatively assess the feasibility of a secondary source of ionization being able to account 
for some of the line ratios. Although no model reproduces the line ratios to within the 
errors, we do find that the composite models better reproduce the run of line ratios than 
photo-ionization models alone. Rand (1998) finds a similar conclusion for NGC 891, where 
composite modeling reproduces the run of [O III]/H/3 vs. [N II]/Hq; and [S II]/Hq;, yet is 
unable to duphcate the [S II] /[N II] ratio at low- 2;. NGC 4302 is not considered due to the 
lack of critical [O III] A5007 data. 
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3.3.1. NGC 5775 

For the composite modeling, we consider separately the cases of filamentary and 
non-filamentary emission. The NE portion of Slit 1 crosses the prominent NE Ha filament, 
previously identified by CRDW as being associated with a shell-like H I structure and 
possibly being involved in vigorous disk-halo interactions. In addition, the SW portion of 
Slit 2 crosses a bright plume of Ha emission immediately adjacent to the very prominent 
SW Ha filament, also a region where intense disk-halo activity is suspected. Since these 
two regions appear associated with current cycling of material from disk to halo, it follows 
that shock ionization may play a greater role in energizing DIG in these areas. The SW 
portion of Slit 1 and the NE portion of Slit 2 each cross regions with relatively less extended 
DIG emission, where disk-halo interactions appear far less vigorous than regions with more 
filamentary features. 

The runs of [O III]/Ha, [S II]/Ha, and [O I]/Ha vs. [N II]/Ha are shown in Figures 
8, 9, and 10, respectively. Figures 8 and 9 also present filamentary and non-filamentary 
regions separately. The runs of these ratios in the matter-bounded S94 models, with various 
values of U labeled, are plotted along with predicted line ratios from the SM79 models 
for Vs — 100 km s~^, Vg = 110 km s~^, and Vg — 100 km s"^ with depletions. We find 
the Vs — 100 km s~^ model is best able to reproduce the observed run of line ratios as a 
secondary source to photoionization. The Vs — 110 km s~^ model has difficulty reproducing 
the run of [S Il]/Ha vs. [N Il]/Ha, and the depleted Vs = 100 km s~^ model, though able to 
duplicate the run of [S Il]/Ha vs. [N Il]/Ha, is unable to reproduce the run of [O IIl]/Ha 
vs. [N II] /Ha at high-2;. 

We find that the run of [O III] /Ha vs. [N II] /Ha seen in Figure 8a for the filamentary 
regions is distinctly different from the trend seen in Figure 8b for the non-filamentary 
regions. Glearly, for a given value of [N II] /Ha, the ratio [O III] /Ha is much greater in DIG 
with a more filamentary morphology. The non-filamentary regions show a run that is very 
similar to the run observed in NGG 891 (Rand 1998). The run of [S II]/Ha vs. [N II]/Ha 
on the other hand, does not seem to depend strongly on the presence of filamentary 
DIG structure. Non- filamentary regions do seem to have greater values of [N II] /Ha and 
[S II] /Ha, in keeping with these ratios' anti-correlation with relative Ha surface brightness 
(see also Figure 12). All data for the ratio [O I]/Ha, being detected only for Slit 2 to 2; = 2.5 
kpc, are plotted in Figure 10. The lack of [O I] /Ha data for higher values of [N II] /Ha 
makes a comparison between various DIG morphologies impossible. 

Gomposite models are plotted in Figures 8, 9, and 10 as dashed lines joining open 
hexagons, which mark various values of log U. For the filamentary regions, 20% of the Ha 
emission at log U — —4.0 arises from shock ionization, increasing to 85% at log U — —5.0. 
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This does not mean the majority of extraplanar DIG is shock ionized. It does indicate 
however that at very high-z (log U = —5.0 corresponds to a height of about z = 6 — 8 
kpc for the filamentary regions and z = 3 — 5 kpc for the non-filamentary), it is possible 
that most of the emission within the filaments arises from shock ionized gas. In the 
non-filamentary case, shocks contribute 16% of the Ha emission at log U = —4.0, rising to 
27% at log U = -5.0. For the [O 1]/Ha vs. [N II]/Ha data of Figure 10, the filamentary 
composite model is plotted over all data points for both filamentary and non-filamentary 
regions. This model and the non-filamentary model are somewhat successful in explaining 
the [O 1]/Hq; vs. [N II]/Hq; run at low-z, in the sense of matching the excess of [O 1]/Hq; for 
a given [N II] /Ho; compared to the S94 model alone. 

The composite models both match the runs of [O III] /Ho; vs. [N II] /Ha. The main 
deficiency in the models is their inability to reproduce the run of [S II] /[N II] at low- 2;. 
The models do not predict enough [S II] emission for a given value of [N II]. The models 
perform somewhat better at high-z, though the spread in data values at high- 2; make a 
trend somewhat difficult to discern. Again, though models with depleted abundances are 
better able to produce an excess of [S II] emission relative to [N II], they do not produce 
enough [O III] emission. 

In any case, it is clear that the filamentary regions require a greater contribution 
from shocks to explain the observed line ratios. The correlation of these features with H I 
shells and enhanced 20-cm radio continuum emission is also suggestive of a larger percent 
contribution to the observed emission in these regions from shock ionization. Enhanced 
[O III] emission can also be produced by other sources of ionization such as TMLs. Rand 
(1998) found that composite photo-ionization/TML models were just as successful in 
reproducing the runs of observed line ratios in NGC 891 as composite models featuring 
shocks. TMLs are expected to occur at interfaces of gas at different temperatures in the 
ISM, such as in shell walls. It is likely then for NGC 5775, where H I shells are associated 
with the more prominent Ha filaments, that TMLs could also account for the enhanced line 
ratios. We have performed cursory modeling of the line ratios for NGC 5775 with composite 
photo-ionization/TML models. Using the TML model of Slavin, Shull, & Begelman (1993) 
with depleted abundances in the mixing layers of hot and warm gas, a hot gas mixing speed 
of Vt = 25 km s^^, and a mixed gas temperature of log T = 5.3, we adopt the same approach 
as the composite photo- ionization/shock modeling. We find that the photo-ionization/TML 
models work nearly as well as the photo- ionization/shock models, though they are unable 
to reproduce the run of ratios for the highest values of [O III] /Ha. 
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3.3.2. UGC 10288 

A preliminary analysis of the diagnostic diagrams for UGC 10288 has revealed no 
discernible trend in the run of [O III] /Ha vs. [N II] /Ha. To attempt to better reveal any 
possible trend, line ratios on either side of the disk have been averaged for each slit to 
estabhsh a run of each ratio vs. \z\. These averaged runs of [O III]/Ha, [S II]/Ha, and 
[O Ij/Ho; vs. [N Ilj/Ho; are shown in Figure 11. The matter-bounded S94 models, with 
various values of U labeled, are plotted along with predicted line ratios from the SM79 
models for Vg = 100 km s~^, = 110 km s~^, and Vg = 100 km s~^ with depletions. Even 
with the averaging, the [O III] /Ha vs. [N II] /Ha run doesn't show an obvious trend: Slit 1 
seems to show a slight trend of an increase in [O III] /Ha with [N II] /Ha, while Slit 2 data 
is clustered near the photo-ionization model. For this reason, we attempt to fit the run of 
[S II] /[N II] instead of [O III] /Ha vs. [N II] /Ha. We adopt the = 100 km s"^ model with 
depleted abundances as the other two models do not produce enough [S II] emission relative 
to [N II]. 

Composite models are plotted in Figure 11 as dashed lines joining open hexagons, 
which mark various values of log U. We find that the composite models best agree with 
the data when shocks begin contributing to emission at log U — —3.5. For the composite 
model, 20% of the Ha emission at log U — —3.5 arises from a 100 km s~^ shock with 
depleted abundances, rising to 45% at log U — —5.0 (log U — —5.0 corresponds to a height 
of about z 1000 pc). 

The composite model successfully replicates the run of [S II] /[N II] at all but very high-z 
where matter bounded photo-ionization models cannot explain values of [N II] /Ha and 
[S II] /Ha greater than unity. Radiation bounded photo-ionization models feature greater 
values of these line ratios; however, predicted runs of [O III] /Ha cannot be reconciled with 
observations. The run of [O I] /Ha vs. [N II] /Ha is also well explained by the composite 
model. The lack of any obvious trend in the run of [O III] /Ha vs. [N II] /Ha makes it 
difficult to assess whether shocks are a viable contributor to the Ha emission. However, 
the high values of [O III]/Ha in Slit 1 do suggest a large contribution from shocks in those 
regions. It should be noted, though, that the galaxy has a fairly quiescent appearance in 
the Ha image with little extraplanar emission, calling into question the notion that shocks 
are permeating the ISM. It would be useful to obtain more [O III] data for this galaxy to 
better diagnose the possible contributions from secondary ionization sources. Nevertheless, 
a clear departure from a pure photoionization model is indicated by the data. 
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3.4. DIG Temperature as the Main Cause of Line Ratio Variations 

In this section, following the approach of HRT, we attempt to determine whether 
variations in line ratios with z can be explained by a change in gas temperature as one 
moves off the midplane. A temperature increase with z provides a simple explanation of 
why all ratios of forbidden lines to B aimer lines are generally seen to increase with z. In 
fact, HRT find for the Reynolds Layer that the increase in [S II] /Ha and [N II] /Ha, and the 
constancy of [S II]/[N II] with z, up to z = 1.75 kpc, can be accounted for by an increase 
in gas temperature alone, without having to invoke a secondary source of ionization. Such 
a scenario requires additional non-ionizing heating for low-density, high-2; gas, such as 
photo-electric heating from dust grains (Reynolds & Cox 1992) or dissipation of interstellar 
turbulence (Minter & Balser 1997). Using line ratio data versus z for NGC 4302, NGC 
5775, UGC 10288, and NGC 891, we attempt to test this scenario by determining gas 
temperatures as well as ionization fractions of constituent elements. 

For collisionally excited ions, the equation for intensity of emission lines is given by 
Osterbrock (1989), 

/.(photons cm-2 s'^ sr'^) = ^ / n,ne^^^^^^^^^^e-(^-/'=^)rf/, (1) 

47r J -/ 4 ■ (jJi 

where and Eij are the collision strength and energy of the transition, cui is the 

statistical weight of the ground level, T4 is the temperature in units of 10^ K, and f^, is the 
fraction of downward transitions that produce the emission line. 

To simplify matters, we assume no change in physical conditions, namely temperature, 
density, ionization fraction, and gas abundance, along a given line-of-sight through the DIG 
layer. Such uniformity of physical parameters within the ISM is clearly not accurate when 
lines-of-sight cross small scale structures such as H II regions or shocks, where temperature 
or ionization conditions may change dramatically. The assumption is valid, however, if 
these parameters are thought of as averages along a given line of sight instead of describing 
localized conditions within the ISM. 

This assumption, along with the use of the collision strengths provided by AUer (1984), 
allows the intensity of [S II], [N II], and [O III] emission to be represented, in units of 
rayleighs, by: 

/67i6(R) = 2.79 X 10^ (I) ^^Tr'-'''e-'-''/^^EM, (2) 
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7a583(R) = 5.95 X 10^ Q |^Tr°-^^^e--«/^^£;M, (3) 



and 

-^5007 



(R) = 4.84 X 10^ ( ^) ^^IMr-^-^s^-^-sr/n^M^ ^4) 
The intensity of Ha emission is given by, 



where EM = J nidi is the emission measure in units of pc cm~^. 

A further simphfication is to assume N+ and H+ have approximately equal ionization 
fractions. This assumption is justified for two reasons. First, photoionization models of S94 
predict that N^/N^ tracks H^/H^ due to similar first ionization potentials (14.5 eV versus 
13.6 eV) and a weak charge exchange reaction. Second, N is not hkely to ionize higher 
than N+ in a pure photoionization scenario due to a second ionization potential of 29.6 eV 
(HRT). A high ionization fraction of O''""'", however, would imply a significant fraction in the 
N"*""*" state as the second ionization potential of oxygen is 35.1 eV. Thus we will be able to 
check the validity of this assumption with the [O III] /Ha data. Due to this simplification, 
of the line ratios available from this data, [N II] /Ha is the best probe of gas temperature. 

The forbidden lines are all dependent on their abundance relative to atomic hydrogen. 

The issue of abundances is a complicated one, the understanding of which is critical for 
interpreting the line ratios. Current evidence suggests that abundances in halo gas may be 
quite different from gas in the disk. Savage and Sembach (1996) find increasing gas-phase 
abundances in the Milky Way from disk to halo, possibly as a result of partial grain 
destruction due to shocks. However, models by Sembach et. al. (2000; hereafter SHRK) 
predict that the ratios of forbidden lines to Balmer lines tend to rise as abundances fall due 
to decreased cooling efficiency. For this reason, it is unlikely that abundance variations can 
explain the increase in ratios of forbidden lines to Balmer lines. However, observations of 
the [O II] A3727 line, to complement [O III] A5007 and [O I] A6300 measurements, would be 
very useful to completely decouple line ratios (which depend on excitation and temperature) 
from variations in abundance. To simplify matters, we adopt abundances that do not vary 
with z. Using galactic gas-phase abundances of S,N, and O of (S/H)= 1.86 x 10^^ (Anders 
& Grevesse 1989), (N/H)= 7.5 x 10"^ (Meyer, Cardelh, & Sofia 1997), and (0/H)= 3 x 10"^ 
(Cardelli & Meyer 1997), along with the assumption of equal fractions of N+ and H+, we 
arrive at expressions relating line ratios to temperature and ionization fraction: 

M ^ 12.2 rr26g-2.18/r4^ (g) 
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- 14 S f 7^0.307 -2.14/r4 (j\ 



and 



_ 1 9 I _SV^ ) Tn-0.119_0.04/T4 

[Nil]" Vh+ZhJ^ 



I H+/H ) 



(8) 



(9) 



The ionization fractions of oxygen and hydrogen are coupled through a charge exchange 
reaction that allows the ionization state of Hydrogen to be determined from the [O I] /Ha 
ratio. We use the ratio of emissivities of the [O I] A6300 line to the Ha line provided by 
Reynolds et. al. (1998a) to obtain, 



[O I] 



1 - 



El 

H 



(^) 



7^1.85 

li g-2.284/T4 QQX 



Gas temperatures are derived from the [N II] /Ha ratio. Assuming this ratio accurately 
reflects temperatures within the medium, line ratios involving the S and O lines are then 
calculated for a variety of ionization fractions. Wc can then sec whether all the ionization 
fractions are reasonable and consistent with our initial assumptions. Plots of temperature 
versus z and calculated ratios overlayed on measured ratios are shown in Figures 3, 4, 
12, and 13 for NGC 4302, UGC 10288, NGC 5775, and NGC 891, respectively. Resulting 
ionization fractions from the temperature modeling are shown in Table 3. 

To first order, a rise in DIG temperature can account for the observed rises in line 
intensities relative to Ha emission, as well as the near constancy in the [S II] /[N II] ratio. 
We find, except in the case of UGC 10288, the ionization fraction of S+ relative to H+ 
remains very close to 0.5 independent of z. This is similar to the result obtained by HRT 
for the WIM. Photo- ionization models performed by SHRK predict a comparable S"*" 
ionization fraction for low-density H II regions. UGC 10288 on the other hand, shows a 
considerably higher fraction of S"*" relative to H"*" as well as some evidence for an increase 
with z. However, the ionization fraction of H"*", determined from points relatively close to 
the midplane, is significantly lower than in the case of NGC 5775 or NGC 891. As a result, 
the range in ionization fraction of 8+ in UGC 10288 is nearly equal to the range in values 
for the other galaxies. This does not explain the observed rise with z, however, and the 
lack of [O I] data at high-z makes it impossible to determine whether the ionization state of 
sulfur changes with z. 
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The [O III] /Ha data, in addition to revealing information about the ionization state 
of O, is also an excellent test of the validity of some of the assumptions used for this 
temperature modeling. The data for NGC 5775 and UGC 10288 suggest that a significant 
fraction of oxygen is in the O^"'' state; at least 15% and in some lines-of-sight nearly 100% 
relative to H"*". The second ionization potential of oxygen is 35.1 eV, while that of nitrogen 
is 29.6 eV. Thus one would expect the ionization fractions of the doubly-ionized species of 
these two elements to roughly track one another in a warm medium. The determined values 
of ionization fraction of O''"''" for these two galaxies imply that a significant percentage 
of nitrogen atoms are in the N++ state, particularly at high-z for NGC 5775. The fikely 
presence of significant amounts of N++ as well as the possibihty of its fraction changing 
with z complicates the use of [N II] /Ha as a tempcratiue indicator. In any event, the high 
values of O^^ ionization fraction suggest that a temperature increase alone is not able 
to fully explain the runs of line ratios observed in NGC 5775 and UGC 10288 since pure 
photo-ionization models predict very low 0^+ fractions. One possible explanation for the 
large ionization fractions of doubly-ionized species is that a secondary ionizing source with 
emission characteristics similar to shocks or TMLs contributes to the layer's ionization. 
Alternatively, the [O III] emitting gas could represent a component that is at a significantly 
higher temperature than the component responsible for the [N II] emission. 

NGC 891 on the other hand, shows considerably lower values for the ionization fraction 
of relative to H+: 13% on either side of the disk and just above 30% in the midplane. 
The sht runs near a bright filament that is not as prominent as those observed in NGC 
5775. These fractions are slightly larger than those predicted by the SHRK model for 
low-density H II regions, though still more than an order of magnitude higher than that of 
the diffuse gas model of SHRK. Thus it seems the run of line-ratios in NGC 891 is better 
explained by an increase in gas temperature with z than for NGC 5775. A slightly warmer 
[O III] emitting component than indicated by the [N II] /Ha ratio could explain measured 
values of [O III] /Ha while keeping the ionization fraction of O"'""'" at levels predicted by the 
models. If this is the case, a secondary ionizing source would not be necessary to explain 
the runs of fine ratios in NGC 891. 

The run of [Sll]/Ha in NGC 4302 is well fit by the increase in temperature reflected 
by the [Nil] /Ha ratio. It would be interesting to obtain [Olll]A5007 data for this galaxy to 
better assess whether line ratio variations can be explained by DIG temperature variations. 

Other factors may influence the line ratios such as afore-mentioned abundance gradients 
(Savage & Sembach 1996) as well as scattered light due to extraplanar dust. The effects of 
dust absorption on extraplanar emission have been shown to be significant for a number of 
DIG halos (Howk & Savage 1999). Modeling of the Milky Way (Wood & Reynolds 1999) 
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and NGC 891 (Ferrara et. al. 1996) suggest that scattered light makes a minor contribution 
to diffTisc halo Ha emission that decreases with z. In NGC 891, the contribution to DIG 
emission from scattered light falls to 10% ai z = 600 pc (Ferrara et. al. 1996). In any case, 
dust-scattered light may be an important consideration for the interpretation of these line 
ratios. 



4. Conclusions 

Deep spectroscopy of four edge-on spiral galaxies has allowed physical conditions 
within their DIG halos to be probed, and has provided information on possible sources of 
ionization of the gas. This study has revealed the following conclusions: 

1. The general trend for the runs of hue ratios for NGC 4302 and UGC 10288 is 
very similar to that of NGC 5775 (Rand 2000) and NGC 891 (Rand 1998): an increase 
in [N II] /Ha and [S II] /Ha with [S II] /[N II] remaining nearly constant with though 
one slit position for UGC 10288 does show evidence for an increase in [S II]/[N II] with z. 
Determinations of halo emission scale height have revealed that the more prominent DIG 
halos also have greater scale height. This is consistent with models of dynamic halos where 
greater levels of star formation in the disk can push material farther off the midplane. 

2. Pure photo-ionization models do not explain the runs of line ratios in NGC 5775 and 
UGC 10288. Composite photo-ionization/shock models can replicate the run of [O III]/Ha 
vs. [N ll]/Ha in NGC 5775, though the run of [S II]/Ha vs. [N II]/Ha is still problematic, 
especially at low- 2;. In addition, we find that the composite models require a greater percent 
contribution from shocks in filamentary than non-filamentary regions to match the run of 
[O IIl]/Ha. A lack of a discernible trend for the [O IIl]/Ha vs. [N ll]/Ha run in UGC 
10288 makes it difficult to determine whether shocks contribute to line emission, though 
the high values of [O III] /Ha in Slit 1 indicate that a shock contribution may be necessary. 
A composite model, though, can account for the run of [S II] /Ha vs. [N II] /Ha. However, 
if it follows the trend seen in NGC 5775 of more disturbed regions requiring a larger shock 
contribution, then the lack of a bright, filamentary halo in UGC 10288 may imply that 
shocks do not play as significant a role in layer's ionization. 

The [O III] /Ha vs. [N II] /Ha run seems to be a key discriminant for shocks. 
Fabry-Perot observations of these 3 lines for brighter halos, such as NGC 5775, are now 
feasible and would provide much more information on the importance of secondary ionizing 
sources with local morphology and height. It is important to emphasize that shocks are not 
the only possible secondary source. Other sources such as TMLs, X-rays, and ionization by 
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cosmic rays may contribute as well. 

3. Another scenario that we consider is a pure photoionization model with an 
increasing temperature with z. By making a few reasonably valid assumptions, the run 
of gas temperature vs. z has been determined in each of the four galaxies. We find that 
an increase in temperature with z can explain the general trend of an increase in the 
ratios of forbidden lines to Balmer lines with including an increase in [O III] /Ha and 
the constancy of [S II] /[N II]. By establishing the run of gas temperature, constraints 
on ionization fractions of Sulfur, Oxygen, and Hydrogen could then be established. The 
[O III] /Ha ratios in NGC 5775 and UGC 10288 generally indicate an exceptionally high 
ionization fraction of 0"'"+, that invalidates the assumption that and H+ have equal 
ionization fractions, indicating that the line ratios cannot be explained by a temperature 
increase alone. The runs of line ratios in NGC 891, where the ionization fraction of 
is still somewhat high, and NGC 4302 are better explained by an increase in temperature 
with ^, though the absence of [O lll]A5007 data for NGC 4302 makes a proper assessment 
difficult. 

Finally, [O l]A6300, [O ll]A3727, and [O 1II]A5007 data for a vertically oriented slit 
would provide additional information on ionization and excitation conditions in these DIG 
halos. 

We thank the KPNO staff for their help in obtaining the data. We also thank L. M. 
Haffner, R. Reynolds, and R. Benjamin for comments in the preparation of this paper. This 
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Fig. 2. — The slit positions for the observations of UGC 10288 overlayed on the Ha image 
of Rand (1996). The image does not show the entire 5' sht length. 
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Fig. 3.— Plots of [NII]/Ha, [SII]/Ha, and [SII]/[NII] vs. z for NGC 4302. The top plot is 
the z dependence of the gas temperature determined from the [N1I]/H« ratio (see section 
3.4). See caption for Figure 12 for an explanation of symbols, curves, and labels used in the 
plots. 
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Fig. 4.— Plots of [NII]/Ha, [SII]/Ha, [SII]/[NII], [OIII]/Ha, and [OI]/Ha vs. z for both 
slit positions for UGC 10288. The top plot is the z dependence of the gas temperature 
determined from the [Nil] /Ha ratio (see section 3.4). See caption for Figure 12 for an 
explanation of symbols, curves, and labels used in the plots. 
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Fig. 5. — Plots of logarithmic Ha intensity vs. z {points) for NGC 5775, with exponential 
least-square fits to the data {dashed lines) for (a) sht position 1 with best fit exponentials of 
Hem = 2300 pc on each side, and (6) sht position 2 with best fit exponentials of — 2500 
pc on the NE side, and Hem = 1800 pc on the SW side. 
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Fig. 6. — Plot of logarithmic Ha intensity vs. z {points) for NGC 4302, with exponential 
least-square fits to the data {dashed lines). The best fit exponentials to the east and west 
sides are Hem — 700 pc and Hem — 550 pc, respectively. 
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Fig. 7. — Plots of logarithmic Ha intensity vs. z (points) for UGC 10288, with exponential 
least-square fits to the data {dashed lines) for (a) slit position 1 with a best fit exponential 
of Hem = 320 pc, and (6) sht position 2 with a best fit exponential of Hem = 460 pc. In 
both cases only emission on the south side of the disk is well fit by an exponential. 
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Fig. 8. — Diagnostic diagrams of tlie line ratios [OIII]/Hq; vs. [Nil] /Ha for (a) filamentary 
regions and (6) non-filamentary regions in NGC 5775. The solid curve represents predictions 
from the matter-bounded photo-ionization model of S94, with various values of ionization 
parameter, log U, indicated by circles. The triangle, square, and pentagon indicate line 
ratios from SM79 shock models with shock speeds of 100 km/s and 110 km/s with standard 
abundances, and 100 km/s with depleted abundances, respectively. The composite photo- 
ionization/shock model is represented by open hexagons connected by dashed lines. The 
hexagons represent ionization parameters for the photo-ionized component of log U = 
—4.0, —4.2, —4.5, —5.0. The fraction of emission from the 100 km/s shock ranges from 20% 
at log U = —4.0 to 85% at log U = —5.0 for the filamentary regions and 16% at log U = —4.0 
to 27% at log U — —5.0 for the non-filamentary regions. 
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Fig. 9. — Diagnostic diagrams of the line ratios [SIIj/Ha vs. [NII]/Hq; for (a) filamentary 
regions and (b) non- filamentary regions in NGC 5775. See the caption for Figure 8 for an 
explanation of symbols used in the diagram. 



-34- 




-35- 




Fig. 10.— Diagnostic diagram of the line ratios [01] /Ha vs. [Nil] /Ha in NGC 5775 for Slit 
2. Symbols are the same as those described in the caption for Figure 8, except we only plot 
the composite model at lowest-^; (the open hexagon is for log U = —4.0). The composite 
model shown is the model appropriate for filamentary regions. 
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Fig. 11.— Diagnostic diagrams of line ratios for UGC 10288: (a) [OIII]/Hq; vs. [Nil] /Ha, 
(6) [SII]/Hq; vs. [NII]/Hq;, and (c) [OI]/Ha vs. [NII]/Hq;. Line ratios on either side of the 
disk have been averaged for each slit. See Figure 8 for an explanation of symbols used in 
the diagram. Open hexagons describing the run of the composite model represent values 
of ionization parameter, log U = —3.5, —4.0, —4.2, —4.5, —5.0, except for the [01]/Hq; vs. 
[NII]/Hq; composite model which shows only the run at low- 2; (log U = —3.5, —4.0. — 4.2). 
The fraction of emission from the depleted-abundance 100 km/s model ranges from 20% at 
log U = -3.5 to to 45% at log U = -5.0. 
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Fig. 12.— Plots of [NII]/Ha, [SII]/Ha, [SII]/[NII], [OIII]/Ha, and [OI]/Ha vs. z for both slit 
positions for NGC 5775. The top plot is the z dependence of the gas temperature determined 
from the [Nil] /Ha ratio. Dashed curves in the line ratio plots are ratios calculated from 
the gas temperature for various ionization fractions. The relevant ionization fractions are: 

(h^/h) [Snj/Ho; and [SII]/[NII], ( ^V/h ) for [OIII]/Hq;, and (^Sl) for [OI]/Hq;. 
Numbers adjacent to the dashed curves in the line ratio plots indicate the value of the 
relevant ionization fraction. 
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Fig. 13.— Plots of [NII]/Ha, [SII]/Ha, [SII]/[NII], [OIII]/Ha, and [OI]/Ha vs. z for NGC 
891. The top plot is the z dependence of the gas temperature determined from the [Nil] /Ha 
ratio. See caption for Figure 12 for an explanation of symbols, curves, and labels used in 
the plots. 
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Table 1. SUMMARY OF OBSERVATIONS 



Galaxy 


Galaxy Center*^ 


Adopted 


Offset along 


Hours of 


Noise*^ 


and Slit 


R.A. 


Dec. 


Distance'^ 


Major axis 


Integration 






(J2000) 


(J2000) 


(Mpc) 








NGC 5775 (Slit 1) 


14'*53"'57!6 


03°32'40" 


24.8 


32" NW 


6.0 


2.2 


NGC 5775 (Slit 2) 








20" SE 


5.0 


2.4 


NGC 4302 


12 21 42.5 


14 36 05 


16.8 


24" S 


4.0 


9.7 


UGC 10288 (Slit 1) 


16 14 25.1 


-00 12 27 


31.5 


30" E 


6.17 


2.5 


UGC 10288 (Slit 2) 








35" W 


5.33 


3.6 


NGC 891 


02 22 33.1 


42 20 48 


9.6 


100" N 


1.5 


16.1 



^The H I kinematic center of NGC 5775 was determined by Irwin (1994). Central positions of 
NGC 4302, UGC 10288, and NGC 891 are from de Vaucouleurs et al. (1991). 

^Adopted distance of NGC 5775 is from Irwin (1994). Distances for NGC 4302, UGC 10288, 
and NGC 891 are from TuUy (1988). 

'^Before spatial averaging. Noise units are 10~^^ erg cm~^ s~^ arcsec"^. 
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Table 2. DIG SCALE HEIGHTS 





Profile 


Hem 


Galaxy 


Description 


(pc) 


NGC 5775 


Slit 1 (NE) 


2300 


NGC 5775 


Slit 1 (SW) 


2300 


NGC 5775 


Slit 2 (NE) 


2500 


NGC 5775 


Slit 2 (SW) 


1800 


NGC 4302 


East 


700 


NGC 4302 


West 


550 


UGC 10288 


Slit 1 (S) 


320 


UGC 10288 


Slit 2 (S) 


460 



Table 3. CONSTRAINTS ON IONIZATION FRACTIONS* 



Galaxy Gas Temperature (^ H"''/h ) ("fFTH") (^fT") 

and Slit (10^ K) 



NGC 5775 (Slit 1) 


0.64- 


-0.95 


0.30- 


-0.60 


0.15-0.70'^ 


n/a 


NGC 5775 (Slit 2) 


0.62- 


-0.95 


0.35- 


-0.80 


0.30-0.75^ 


0.65-0.85 


UGC 10288 (Slit 1) 


0.60- 


-1.05 


0.55- 


-1.00 


0.50-0.80*^ 


0.45-0.65 


UGC 10288 (Slit 2) 


0.60- 


-0.95 


0.50- 


-1.00 


0.50-0.80^ 


0.30-0.60 


NGC 4302 


0.66- 


-1.00 


0.35- 


-0.60 


n/a 


n/a 


NGC 891 


0.65- 


-1.05 


0.40- 


-0.65 


0.13-0.30 


0.55-0.80 



^Values of ionization fractions are valid only for the case of a pure 
photoionization model with an increasing temperature with z. See text for other 
assumptions used to determine these values. 

^Nearly all data points fall within this range for NGC 5775 and UGC 10288. 
However, some measured values indicate an ionization fraction of O"'""'" relative 
to H+ that is nearly 1.0. 



